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ABSTRACT 

Measurements of distributions associated with the pair production of top 
quarks at the LHC can be used to constrain (or observe) the anomalous chromo- 
magnetic dipole moment (k) of the top. For example, using either the tt invariant 
mass or the pt distribution of top we find that sensitivities to \k\ of order 0.05 are 
obtainable with 100 fb~ l of integrated luminosity. This is similar in magnitude 
to what can be obtained at a 500 GeV NLC with an integrated luminosity of 50 
fb~ l through an examination of the e + e~ — ► tig process. 
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1 Introduction 



The Standard Model(SM) has provided a remarkably successful description of almost all 
available data involving the strong and electroweak interactions. In particular, the discovery 
of the top quark at the Tevatron with a massjl|, m t = 175 ±6 GeV, close to that anticipated 
by fits to precision electroweak data[|| is indeed a great triumph. However, we know that 
new physics beyond the SM must exist for many reasons particularly those associated with 
the fermion mass generating process. Since the top is the most massive fermion, it is believed 
by many that the detailed physics of the top quark may be significantly different than what 
is predicted by the SM. In this scenario, the top provides a window into the new physics 
which lies beyond the electroweak scale. This suggestion makes precision measurements of 
all of the top quark's properties absolutely mandatory and will require the existence of top 
quark factories. 

One of the most obvious and easily imagined scenarios is one in which the top's couplings 
to the SM gauge bosons, i.e., the W, Z, 7, and g, are altered. In the case of the electroweak 
interactions involved in top pair production in e + e~ collisions, the lowest dimensional gauge- 
invariant operators representing new physics that we can introduce take the form of dipole 
moment-type couplings to the 7 and Z. In the case of strong interactions, the subject of 
the present work, the corresponding lowest dimensional operator conserving CP that we 
can introduce is the anomalous chromomagnetic moment «|3|, M. On the otherhand, the 
corresponding chromoelectric moment, k, violates CP. In this modified version of QCD for 
the top quark the tig interaction Lagrangian takes the form 



where g s is the strong coupling constant, m t is the top quark mass, T a are the color generators, 
is the gluon field and q is the outgoing gluon momentum. Due to the non-Abelian nature 
of QCD, a corresponding four-point tigg interaction, proportional to k and/or k, is by 
necessity also generated. 

Perhaps the most obvious place to probe for anomalous top couplings is at hadron col- 
liders. It is clear that the existence of a non-zero value for n (and/or k) would lead to a 
modification in both the gg — > ti and qq — > ti subprocess cross sections at these machines. 
The general expressions for these parton level cross sections are given in Atwood et al. ||. 
Here we note only that the qq subprocess has a quadratic k dependence while that for the 
corresponding gg subprocess has a quartic dependence on k. In our discussion of anomalous 
top couplings at the LHC, we will ignore for brevity the possibility of a non-zero k. Obvi- 
ously, the observation of the CP- violation induced by non-zero k is a more sensitive probe 
for the anomalous chromoelectric moment of the top than the kinematic distributions we 
consider below. 
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Figure 1: Cross section for ti production as a function of n at the Tevatron for m t = 175 GeV. 
The dotted(dash-dotted) curve is the qq(gg) contribution and the solid line is their sum. MRSA' 
parton densities were assumed. The horizontal dashed bands correspond to the la world average 
top pair cross section obtained by CDF and DO. 
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2 Effects of Anomalous Couplings 



At the Tevatron, it has been shown|| that for small values of \k\ < 0.25, a range consistent 
with the current total cross section measurements [[!]] by both CDF and DO, the dominant 
effect of anomalous chromomagnetic moment couplings is to modify the total cross section for 
top pair production with little influence on the shape of the various distributions. Figure |l| 
compares the K-dependent cross section with the world average of that obtained by the CDF 
and DO Collaborations. 

The essential reason why the various top quark kinematical distributions are not much 
influenced is that top pair production at the Tevatron is dominated by the invariant mass 
region near threshold. Since, as is well known, the effects of anomalous couplings grow with 
the parton center of mass energy one sees little influence at these energies. The significantly 
larger partonic center of mass energies accessible at the LHC allows us to probe beyond this 
threshold region so that much higher sensitivities to a possible non-zero k can be obtained. 
This is particularly true for the various kinematic distributions. 

As a result of the subprocess dependencies on k it is clear than any of the (unnormalized!) 
differential distributions for a generic observable, O, can then be written as 

% = i: K n 9n (o) (2) 

where g n {0) are a set of calculable functions which have been completely determined to 
lowest order in QCD by Atwood et a/. ||. The QCD/SM result is just the familiar term with 
n = 0. Of course, the total cross section is also a quartic polynomial in k. The behaviour of 
the two individual contributing subprocess as well as the total cross sections under variations 
of k at the LHC are shown in Fig. |2|. Unlike the Tevatron, the gg initial state dominates the 
top pair production cross section at the LHC. A reasonable sensitivity to k is again observed 
in the total cross section as it was for the Tevatron. However, as discussed in Ref. 0, unless 
the theoretical and systematic uncertainties are well under control, a measurement of o~ t t at 
the LHC will never do much better than to constrain \k\ < 0.10 — 0.15. To further improve 
on this limit we must turn to the various top quark kinematical distributions. 



3 Analysis 

As has been shown elsewhere 0], the pt and pair invariant mass (M u ) distributions for top 
quark pair production at the LHC are highly sensitive to non-zero values of k. Figures ^a and 
|3]c show the modifications in the SM expectations for both da/dM tt and da/dp t , respectively, 
for different values of k. Perhaps more revealingly, Figures |3|b and ^d show the ratio of the 
modified distributions to the corresponding SM ones. We see the important results that a 
non-zero k leads to (i) enhanced cross sections at large pt and M tt and (ii) the shapes of 
the distributions are altered, i.e., the effect is not just an overall change in normalization. 
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Figure 2: Cross section for ti production as a function of k at the LHC for m t = 180 GeV. The 
dotted(dash-dotted) curve is the qq(gg) contribution and the solid line is their sum. MRSA' parton 
densities were assumed. 



5 




Figure 3: (a) ti invariant mass distribution at the LHC for various values of n assuming rrit = 180 
GeV. (b) The same distribution scaled to the SM result, (c) ti pt distribution at the LHC and 
(d) the same distribution scaled to the SM. In all cases, the SM is represented by the solid curve 
whereas the upper(lower) pairs of dotted(dashed, dash-dotted) curves corresponds to k =0.5(-0.5), 
0.25(-0.25), and 0.125( -0.125), respectively. 
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This is contrary to what was observed in the Tevatron case where both da/dM u and da/dp t 
were essentially just rescaled by the ratio of the total cross sections. Clearly, data on these 
two distributions at the LHC can lead to significant constraints on k or observe a non-zero 
effect if k is sufficiently large. In Ref. ||, the cosO* and rapidity distributions were also 
examined but they were found to be less sensitive to non-zero k than the dramatic effects 
shown in Fig. |3|. 
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Figure 4: Sample histograms of top quark data generated for the LHC, assuming 100 fb 
of integrated luminosity. On the left (right) is the top pair invariant mass (p t ) distribution. 
MRS A' parton densities and m t = 175 GeV have been assumed. 

How sensitive are these distributions to non-zero k and what bounds can be obtained at 
the LHC? In order to answer these questions, we follow a Monte Carlo approach. We begin 
by generating 100 fb~ l 'data' samples for both distributions assuming the SM is correct. To 
be specific, since the next to leading order(NLO) expressions for these distributions in the 
presence of anomalous couplings do not yet exist, we use the leading order results rescaled 
by the NLO/LO cross section ratios for both subprocesses as effective i^-factors to obtain a 
rough estimate of these higher order effects. Sample histograms of this appropriately rescaled 
'data' are shown in Fig. f|. Note that there are 37 bins in M tt and 22 bins in p t of varying 
sizes essentially covering the entire kinematically allowed ranges. Bin sizes are adjusted to 
partially conform to changes in resolution and declining statistics as we go to larger values of 
either kinematic variable. In addition to the usual statistical errors, we attempted to include 
some estimate of the systematic point-to-point errors. These were added in quadrature to 
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the statistical errors. Thus, neglecting the overall normalization uncertainties, the error in 
the number of events(iVj) in a given bin(i) was assumed to be given by 

SN^iNt + aN?} 1 / 2 (3) 

with the parameter a setting the a priori unknown size of the systematic error. Note that we 
have made the simplifying assumption that the magnitude of a is bin independent. The total 
error is thus generally systematics dominated. With these errors the Monte Carlo generated 
data was then fit to the known functional form of the relevant distribution: 

% = fj:^9n{0) (4) 

au n =0 

where / allows the overall normalization to float in the fit and the g n were those appropriate 
to either the p t or M tt distributions. 

The results of this analysis are thus a set of 95% CL allowed regions in the f — K plane 
for various assumed values of the anticipated size of the systematic errors. These can be 
seen in Figure [5[ Here we see that for systematic errors of reasonable magnitude the value 
of k is constrained to lie in the range —0.09 < k < 0.10 from the M a distribution and 
—0.06 < k < 0.06 from the corresponding p t distribution. Note that the correlation between 
/ and k is much stronger in the case of the M tt distribution. Increasing the integrated 
luminosity by a factor of two will not greatly affect our results since the errors are systematics 
dominated. Combining the results of multiple distributions in a global fit to k will most likely 
result in even strong bounds. 

Using Figure |6| we can make a direct comparison of the bounds obtainable on k at the 
NLC by using the process e + e~ — > tig as discussed in Ref.Q with those from the LHC 
analysis above. In these Figures the influence of R is also shown. These NLC results were 
obtained by fitting the spectrum of very high energy gluon jets produced in association 
with top pairs (above some cut, E™ m , used to avoid contamination from the radiation off 
final state 6-quarks in top decay). Only statistical errors were included in the analysis. 
The resulting bounds are essentially statistics limited. We see from these Figures that the 
constraints on k from the \/s=500 GeV NLC with an integrated luminosity of 50 fb^ 1 are 
only slightly better than what is achievable at the LHC from the top pair's p t distribution. 
The constraints tighten at the 1 TeV NLC. Clearly the LHC and NLC have comparable 
sensitivities to the anomalous chromo magnetic moment of the top. 
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Figure 5: 95% CL two parameter (/, k) fits to the invariant mass(left) and p t (right) dis- 
tributions at the LHC for a 175 GeV top quark assuming the MRS A' parton densities for 
different assumed values of the systematic errors parameterized by a. From inside out the 
curves correspond to a— 0.03, 0.05, 0.10, 0.15, 0.20, and 0.30, respectively. 
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Figure 6: 95% CL allowed regions in the k—k plane obtained from fitting the gluon spectrum. 
On the left the fit is for gluon jets above E™ m =25 GeV at a 500 GeV NLC assuming an 
integrated luminosity of 50(solid) or 100(dotted) fb~ x . On the right is the case of a 1 TeV 
collider with £™ m =50 GeV and luminosities of 100(solid) and 200(dotted) fb~ l . Note that 
the allowed region has been significantly compressed downward in comparison to the 500 
GeV case. 
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